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ABSTRACT 


The major environmental Ultra Low Frequency (ULF), magnetic 


_noise sources observed on the airborne AN/ASQ-81(V) optically 


pumped magnetometer are electromagnetic micropulsatrons, geo- 
logical structures and ocean wave noise. 

Experimental determinations of the amount of magnetic noise 
produced by ocean surface gravity waves were made. The author's 
current research designed to remove the wave noise and increase 


the probability of Magnetic Anomaly Detection (MAD) is summarized. 
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I. INTRODUCTION 


Magnetometry has been used as a classical geophysical pro- 
Epecting tool for years (Dobrin, 1952). The principleg of 
Magnetic anamoly analysis has also been applied to Anti- 
Submarine Warfare (ASW) with a limited amount of success for 
submarine localization. 

The sensitivity of magnetometers has been continually 
improving so that it is now possible to observe, from an air- 
Beart, fluctuatious of toa nano Tesla/hertz (nT/Hz), imposed 
On the earth's constant field of approximately 50,000 nT. 

Land stationed superconducting magnetic sin@sors have demon- 
strated a sensitivity of ie nT/Hz. The advent of airborne 
superconducting magnetometers would greatly increase detection 
ranges$ however, the extent of improvement depends in detail 
on the amplitude and frequency of the environmental noise which 
is above the noise floor of the instrument. 

Airborne Magnetic Anomaly Detection (MAD) for ASW is done 
in the .04 to 2 Hz frequency band. Ocean surface 
in the .03 to 10 Hz frequency band (Silvester, 1974) are con- 
tributors to the magnetic noise spectrum. The other two major 


environmental noise sources are electromagne 
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geological noise. 
It is generally assumed that noise due to ocean gravity 
Waves is not "loud" enough to be of any consequence. However, 


experimentally it was observed that noise levels induced by 


typical sea condition in the mid latitude oceans during the 


Bereater than the noise floor of airborne magnetometers. The 


TiO Se THEO RY. 


The problem of determining the electric current induced 
Eby deep water gravity waves at the ocean surface and the resul- 
tant periodic fluctuations in the magnetic field has been de- 
scribed by Weaver [1965]. The magnetic field fluctuations 


above and below the ocean were obtained by solving Maxwell's 


Equations for the magnetic flux density vector in air and water 
and applying electromagnetic boundary conditions at the sea 
surface. 

The movement of the sea surface through the earth's mag- 
netic field induces electric eddy currents in the sea water. 
In deep water, assuming linear ocean wave theory, the water 
particle motion is periodic and circular, as is the induced 
Magnetic field. This process becomes more complicated as one 
decomposes the sea surface spectrum into waves of various 
frequencies which may be propagating in various directions 
With respect to the magnetic field vector. Further complica- 
tions arises as higher order wave theories are applied to de- 
scribe the water particle motion and the waves propagate into 


shallow water. 


AS MATHEMATICAL MODEL 
The mathematical model of the induced field was developed 
by Weaver [1965], assuming linear wave theory to describe the 


Water particle motion. 


It was derived for a monochromatic ocean surface gravity 
Wave, propagating in deep water, in the X direction. The 
development is summarized below. 

The coordinate system formed a right hand triad as follows: 

X-axis, was taken in the direction of wave propagation. 

Y-axis, was taken parallel to the wave crest. 

Z-axis, was taken downward into the sea, Z=0 at the 
still water level (SWL). 

It was assumed that the earth's magnetic field was constant 
in the region, and the wave crest was infinitely long. The 
model became two dimensional as all magnetic gradients in Y 
were negligible. 

The earth's vector magnetic field, Bo; was expressed as 


™ 


oe =(B,|(cosT cos @ tc - cos[ sin Of + sw] EE) 


nan an A 


Where i, j, k were unit vectors along the three coordinate 
mors X, y, z. I, was the magnetic dip angle in the vertical 
plane, relative to the horizontal. 8, was the inclination 


in the horizontal plane of the direction of wave propagation, 


Ww 
. 


clockwise from the magnetic meridian (see Figure 1 


The field vectors satisfied Maxwell's equation in the form 


“OR : (1) 
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Volts 


_— 
where E was the electric field vector ( /Meter), B was the 


induced magnetic flux density (Tesla), H, was the permeability 


ef free space (umx107? Henry meter), Oo was the electrical con- 


Siemen 


= 
ductivity ( /Meter), V was the water particle velocity 


bts fee). eo was: the permittivity (PePaes jeter). and 


vector ( 
t was time (seconds). 
The deep water particle velocity, described by linear wave 


theory, was of the form 


( (Wt -mXx) -my 
saa AON oF 
V=2@wCic+k)2 


where a was the wave amplitude (Meters), w=2n/f was the an 


Radians 


frequency of the wave (° /sec), m=2n/L was the wave number 


Radians ; 
( /Meter), and t was time (seconds). 
. => . . . . 
Since V was harmonic in x and t, the magnetic and electrical 


= 
field resulting from V has a similar structure, thus 


— a L(wt-mx) 
b= bey Z- ; (4) 
-~ = L(wt -mX) 
E= Ec) a 7 . (5) 


—_ 
Es 
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The B field and the 


Oba), 
BEG Lt Dy (3) > (6) 


field are non-divergent which implied 
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The differential equation in water satisfied by B was 
obtained by assuming B<<B, and we<<o. Then, taking the curl 
of equation (2) and using the non-divergent portion of the 
vector identity, Vx (TxA) = ATCA + V(vV-A), the following was 


obtained: 


CO 
ww 


~ V*B = OVKE+ 6H X(VX Bz) : 


Further simplification of equation (8) was obtained by the 
elimination of Vx= by equation (1). From equation aay was 
eliminated. Then if x? = LU,@9, the following equation results, 

Ob > 2 

Sz = (% mM )by a ye (CBuwo-cBait+R)g a (9) 
In region Z>0, equation (9) was solved for the Z component 


by assuming 


b 


FA oe. ee 
w= oa" 2 4 a.m (BW iBug le 7 


y) ) 
where Q was an undetermined coefficient. 
im air €>>0, thus equation (2) Became, 
B. ge 
=~ Fe 1 
Vx== ES ; (10) 
Z Jt 


the curl of euation (10), when E was eliminated by equation 


(1), resulted in 


CT es 
pies pt tea) bis) ; 
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Confining the sea surface to scales of motion described by 


gravity waves implied that m>>y eu, thus 


3 


VY be i 
¢ rae (12) 
73 

In the region Z<0 the solution for the Z component of 


equation (12) was of the form 
Nate! in Cis) 
bg) = Ke 


where R was an undetermined coefficient. 
The magnetic field vector was continuous across the air-sea 
interface, thus, the following boundary condition at Z=0, existed: 


Dy az (Z>0) 


b, WaTeR (z =O) 


ob; pin (2 =2) Db; waren (Z20) 


oy O4 (14u.d) 


(l 


Applying the above boundary values to equations (9) and (13) 


yielded 


R=Q+ am(Bw - Rats) ) 


and 


oO 
wm 


2. . 
pike 2 ~O (H+ mit) - am (Boo =e Bow) ; (16 


EBquataons (i5) and (16) lead to 


Q- —Z ma ( Race) -¢ Bey? eas 
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and 


Re am (Bete) -e Bag sent) my 
ig 


Thus the Z component of the induced magnetic field vector 


(18) 


in air was yf 2 
ee 
Da(3) = = Qn m (Bato Beg Ow + pm ) “- ym) 2” 
x 


,) 


< 


In order to further simplify equation (19), =<<l was true 


for gravity waves resulting in 


; asx | | 
bay . Aim (Bw -iBay)a"9 
From equation (6) and (20) it was seen that 


ob psi (21) 
xCy) 36? 
It was convenient to measure the earth's field with a 


total field magnetometer. Thus, 


=2) eb cD COS Tees Etbiyy sin T (22) 


joteere 


When equations (20) and (21) were used to simplify equation 
(22), and only the perturbations were considered, it can be 
shown that the magnitude of the observed magnetic field ata 


Mespecific frequency would be, 


b(f) = Te JB ((coszeced (sw ty) 2F (23) 
| ZO , 
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The deep water wave length (L) of a surface wave is related 


to the wave period (T) by, 
Zz 
ee gl” (24) 
Ait 


where g is the acceleration due to gravity. 


Eliminating m, from equation (23) by the use of equation (24) 


leads to, 


z 
fie /8 SHoTa (2 Mone Ly Pre all 
b(t) = ZEA |Bl(est ese) +{sinT) eo ard . es 


Equation (25) was used to model the magnetic field spectrum 
when the directional ocean wave spectrum was known. 
Using an airborne magnetometer to measure the magnetic 


Spectra, there was a doppler type frequency shift, described as, 


f ' i) (26) 
f= fi - (2fu/g) &s #) 
(see Appendix C). Where f' was the wave frequency observed on 
the aircraft (Gees f was the actual wave frequency, u was 


the velocity of the aircraft Gas 


/sec), and w was the angle 
Between aircraft motion and direction of wave propagation (see 


mesure 1). 
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Pil. THE EXPERIMENT 


Magnetic data was collected by an airborne magnetometer 
flown over the sea surface in order to verify the theory pre- 
viously presented and to further understand the observed environ- 
ment in the 0 to 2Hz band. The AN/ASQ-81 magnetometer, aboard 
the P-3C Orion aircraft, is capable of measuring the signal 
produced by ocean gravity waves which are frequency shifted 
due to relative motion between the wave and the aircraft flying 
at 100 meters altitude. The analog voltage time series output 
©f the magnetometer was recorded on an HP3960, FM tape recorder 
for post flight analysis. 

The fluctuation of the sea surface below the aircraft was 
determined by two techniques. The analog voltage time series 
output of the APN-194 radar altimeter, also aboard the aircraft, 
was recorded. 

Independent of the aircraft data the ocean surface wave 
climate was determined from analysis of the Fleet Numerical 
Weather Central (FNWC), spectral Ocean Wave Model (SOWM) and 
Measured wave spectrum from National Oceanic and Atmospheric 
Administration (NOAA), Environmental Buoy, NR-6003 (EB-17), 
mocated at 52.2N, 155.8W from the airborne data. 

The airborne data was collected in vicinity of 52°N, 156°wW, 


around the buoys' reporting period, 00Z, 10 Feb 78. 


mM. ALTRBORNE DATA COLLECTION 
The AN/ASQ-81 employs the atomic properties of optically 


pumped, metastable helium atoms to detect variations in total 


ony, 


mesure 2. Region of data collection. (X)-FNWC Grid Points, 
(A)NOAA Buoy, (---) First Leg of Aircraft Data, 
(--.--) Second Leg of Aircraft Data. 
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field magnetic amplitudes. The output of the sensing element 
is converted to an analog voltage which is bandpassed and scaled 


~llpy as a function of time. 


to display field variations (Order 10 
Instrument specifications are described in Appendix A. 


: fe) fo) ; 
Perpendicular legs, 270° True and 000° True, were flown at 


approximately 60 meters altitude and 180 knots (91.44 M/s) for 
30 minutes each. (See Figure 2) The magnetic data was collected 
gn the .04 Hz to 2 Hz band, thus, eliminating the earth's con- 
Stant field and any micropulsations above 2 Hz. Aircraft motions 
were compensated for by the magnetometer using baroaltimetry, 
selective filtering and accelerometers. No variations in the 
magnetic field due to the aircraft were expected due to the 

above real time processing and the straight and level flight 

paths flown. 

The instantaneous aircraft altitudes were recorded simul- 
taneously with the magnetic data from the APN-194 radar altimeter, 
which responded linearly to apparent altitude changes due to 
long period surface gravity waves. The wave frequency was 
shifted due to aircraft motion (see equation (26)). Beam 

) characteristics, 60° cone, were such that at 60 meters the 
diameter of the radar spot on the surface was approximately 70 
Meters. In order to resolve wave amplitudes from 60 meter 
altitude, the wave period must be greater than 10 seconds. 
The frequency and amplitude errors attributable to the size 
of the radar spot were considered in the analysis. 

The AN/ASQ-81 and APN-194 signal pick-off points are 
medescribed in Appendix A. 


1c 


PB. SPECTRAL WAVE DATA COLLECTION 


The experiment was conducted around synoptic reporting 


period 00Z, 10 Feb 1978 in vicinity of NOAA Buoy, EB-17, so 

Bithe Spectral information computed by the Buoy for the 00Z time 

period was useful in interpreting the aircraft data. 

In addition to Buoy data the SOWM output was used in order 

) to infer directionality and improve the spectral wave estimate 
Bn the region of the flight. 

The Buoy is designed to report meterological and oceano- 
graphic conditions at the sea surface every 3 hours. The wave 
conditions are measured for 20 minutes by accelerometers and 

reported for the Buoy location as frequency and amplitude 


Mm (Steele and Johnson, 1977). The Buoy data of 00Z, 10 Feb 1978 


is reported in Table I. 


Frequency (Hz) Period (sec) Amplitude (meters) 
06 L6.. 7 wOOie2 
Oy 14.3 27574 
0:8 L236 1.6064 
09 cae ee 1.9604 
sal 1OE0 ds 352 
tod ciel dies 4S 
v2 843 LLL 
213 Verdi +9083 
| 14 Toe {7617 
| sus 6.7 6543 
aialile 6.3 atte? 33 
Sa 69 ~4511 
ws 5446 s4212 


Table I 
NOAA Data Buoy 00Z, 10 Feb 1978 
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The FNWC, 


SOWM uses the analyzed surface pressure field 


to determine wind velocities and then builds, radiates and dis- 


Sipates waves 


fo-3). Model 


and direction 
frequency can 


SOWM frid 


as prescribed by wave theory (FNWC Tech Note No. 
output is energy as a function of frequency 

at each grid point. The wave energy (E) ata 
be converted to wave amplitude (a) by a= E-Af. 


coverage of the data area is shown in Figure 2. 


The model output for the closest grid point is summarized in 


Table II. 
DIRECTION FROM 283° chicas 343° oe? 
COT) 

PERIOD (SEC) AMPLITUDES (METERS) 

m6. 4 Bai 7 0 S78 By, 

iy 93 578 4O8 1.002 £578 

13.89 708 993 1.082 au7 

2.82 nigeiseeds 0 Tees 1.686 

12.04 5227 ieee 10 59 959 

10.99 ely, Le5eu 914 1.192 
8.71 1.293 1.493 944 1.436 
a5 944 i. 356 997 1.27% 
Tas $851 2973 987 914 
eats 409 659 5 Bey; 

Table II 


FNWC, SOWM Data, 03Z 10Feb78. Grid point located at 


52.357N, 


154.099W. Amplitudes (meters) 


2a 


IV. DATA ANALYSIS 


me TIME SERIES ANALYSIS 
The general form of the magnetic and altitude time series 


is shown in Figure 3, which is a sample of each time series. 


The time series were digitized at .1 second intervals, 
resulting in a Nyquist frequency of 5 Hz. Any spectral power 
above 5 Hz would be folded back into the spectral estimate 

fe (alliasing). The magnetic data, due to the low pass, 2 Hz, 
characteristics of the instrument, alliased only instrument 
| noise which was at a low level. It was assumed that no sig- 
nificant alliasing of the altitude data occurred due to the 
altimeters 0.1 second RC constant which results in a 1.59 Hz 
cutoff frequency Co, =) / RC). 

Prior to computing power spectral estimates, the data was 
Calibrated to insure constancy of unit magnitudes. The magnetic 
data voltages were linearly calibrated to nT and the altimetry 
voltages were linearly calibrated to meters after the flight. 

The data was prepared for the Fast Fourier Transform (FFT), 
(Bendat and Piersol, 1971) by forming 16 ensembles, 102.4 seconds 
each. Each ensemble was then linearly detrended. The detrending 
removed significant variance from the frequencies below .001 Hz. 
The cosine function was used to taper 10% of the time series 
by multiplying the first and last 51 points of the time series 
By the cosine. Points 1 and 1024 become zero and points 52 
mmr 973 are as digitized. The tapering procedure is done in 


order to prevent non-real low frequency power leakage, due to 
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Horizontal 


Sample of magnetic time series (top) and 


altitude time series (bottom). 


Pigure 3 . 


Vertical scale is lmm/volt. 


scale is lmm/sec. 
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ensemble averaging, into the spectrum. The variance removed 
in the time series taper is restored to the power spectra after 
fhe FFT (Bendat and Piersol, 1971). 

The FFT computes the real and imaginary components of the 
Fourier Transform. It was then possible to compute the sample 
power spectra of each signal by ensemble averaging. The normal- 
ized cross power spectra, defined as coherence and the phase 
of the cross power spectrum between signals were also computed. 

Each of the two flight courses were treated as separate 
data sets. Note in Figure 4 the power spectra of the altimeter 
(mM? /Hz) and magnetometer (nt? /Hz) and the coherence and phase 
as a function of frequency. The numerical results are summarized 


in Appendix D. 


B. SPECTRAL WAVE ANALYSIS 

The analysis of the FNWC, SOWM directional wave spectra 
and EB-17 ocean wave amplitude spectra in order to describe 
the magnetic field variations observed aboard the aircraft, 
mBequired fitting the large scale information to the microscale. 
This change of scale was feasible due to spatial uniformity 
of the ocean wave spectra over the area and frequency range 
investigated. 

The ocean wave climate in an area can be described by that 
at a point over a period of time. EB-17 reported the waves 
during a 20 minute period to have the spectral shape listed in 
Table I. The FNWC ocean wave amplitude and directional spectra 


mad a similar shape listed in Table II. 
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The directionality of wave components at EB-17 was deter- 
mined by group velocity analysis (Silvester, 1974) which pro- 
pagated the FNWC, SOWM data at nearby grid points to the buoy 
position. The greatest weight was on the grid point closest 
to the position of the buoy. The grid point was located at 
52.357N, 154.099W, 60 miles from the buoy. The analyzed buoy 


directional amplitude spectra is reported in Table III. 


PERIOD (SEC) AMPLITUDE (METERS) DIRECTION FROM (°T) 
1627 #6322 340° 
pies .7574 340° 
12) 25 1.6064 315° 
telee tt 1.9604 315° 
eres 1.9312 320° 

owe 1.5419 395° 
8.3 1.0116 00” 
ee) .9083 320" 
Tal 27919 330° 
67 £6543 330° 
6.3 4 Si3 320° 
5.9 4511 s00° 
5.6 4212 a0 
Table III 


Analyzed directional wave data at NOAA buoy EB-17, 00Z 10 Feb 78. 


It is important to note that the NOAA data and the FNWC 
data have similar spectral shapes although amplitudes are not 
equal, verifying scale analysis assumptions. The buoy reports 
a spectral peak at approximately 11 seconds. FNWC predicts a 


Spectral peak at approximately 12 seconds. 
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The observed ocean wave climate on the aircraft, although 
similar, was not the climate observed by EB-17 or predicted 
by SOWM. Each of the ocean wave data sources had errors asso- 
ciated with them, which would contribute to the differences 
in the data. Scale considerations, length of airborne data 
record, modeled grid point locations, buoy location and flight 
path, further contribute to the uniqueness of the separate 
data sources. 

Due to the multidirectional nature of the SOWM data and 
its proximity to the flight path, it was assumed that the 
closest FNWC data point (Table II) most closely resembled the 
ocean wave climate observed on the aircraft. The SOWM wave 
information was applied to equation (25) resulting in an ex- 
pected magnetic amplitude spectrum. The expected frequency 
shift due to the aircrafts motion relative to the waves was 
computed by applying equation (26). To be consistent with the 
power spectral density computed from the aircraft data the 
amplitude spectrum described above was squared and normalized 
by multiplying by the frequency interval of each frequency 
estimate. The results are summarized in Figure 5 and numeri- 


Cally in Appendix D. 
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V. DISCUSSION 


Baker and Graefe (1968) showed experimentally that the 
theory reviewed in Section II describes the magnetic field 
induced by ocean surface gravity waves. Comparison of the 
aircraft magnetic power spectra, Figure 4, and the expected 
magnetic power spectra computed from SOWM data, Figure 5, 
further verifies the electromagnetic theory described above. 


1gure 4a on 


a | 


Figure 6, which is the superpositioning of 
Figure Sa and Figure 4b on Figure 5b,illustrates that althou 
the aircraft and wave data magnetic spectra are not the same 
they are similar in amplitude and frequency. The differences 
observed in the low frequency end can, in part, be attributed 
to the assumptions necessary to the formulation of each magnetic 
power spectrum and the assumption that each data source is 
stationary in time and space. 

The FNWC data, although close to the flight path, was not 
coincident with the aircraft track (see Figure 2). Furthermore, 
determining the barometeric surface pressure field over data 
Sparse areas, such as the Gulf of Alaska caused inaccuracies 
in the SOWM wave height, frequency and directional information, 
which could have contributed to some of the differences observed. 

When studying the magnetic spectrum computed from the air- 
Graft data it cannot be assumed that the only source of noise 
Was the ocean gravity waves. Although no significant electro- 
Magnetic micropulsations in the 0 to 2 Hz frequency band, were 


expected, they cannot be ignored. 
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Production of magnetic noise by the aircraft's motion and 
onboard equipage was considered. Ideally, most aircraft pro- 
duced noises were compensated for in the magnetometers real 
time processing scheme and the flight was straight and level 
on each leg which would further reduce noise. However, expe- 
rience shows that the noise sources described above are not 
completely removed from the spectral analysis when long data 
mecords are required. 

Note in Figure 6a that both the aircraft and the analyzed 
SOWM magnetic power spectra indicate a broad coincident frequency 
band. In Figure 6b it is seen that a change in the flight path 
causes the band width of coincidence to compress, shift to 
higher frequency, and two distinct maxima are observed in the 
aircraft spectrum. One of which is nearly coincident with the 
SOWM magnetic spectrum. This indicates that a significant mag- 
metic noise source is external to the aircraft and ocean surface 
wave dependent. 

Most significant, is the information implied from the analysis 
©f Figure 4 and Figure 6. It is seen (Figure 4) that the co- 
herence is high and the phase relationship is good between the 
Magnetometer and the altimeter in a small frequency band (.2 to 
-/ Hz). The trend observed in the phase relationship over that 
Band can be attributed to signal delays imposed by the electrical 
@haracteristics of each instrument. 

This highly correlated frequency band falls in the same 


band in which the ocean gravity wave noise is expected. The 


3 


mitimeter specifications indicated that it was accurate, + 1 
foot and has a linear response to about 1.5 Hz. Thus it can 
be concluded that the altimeter senses the sea surface flucuations. 
Note the large spectral peak below .1 Hentz on the alti- 
meter power spectra of Figure 4. The peak a approximately 
one order of magnitude greater in the data of the first flight 
leg (Figure 4a) when compared to the second (Figure 4b). This 
peak was due to an intermittent square wave caused by malfunc- 


tion of the altimeter. The square wave does not significantly 


degrade the analysis presented above. 
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VI. APPLICATION 


It was shown in section V that the altimeter can be used 
to describe the instantaneous fluctuations of the sea surface 
below the aircraft which are directly related to a portion of 
the magnetic noise. Thus, the high correlation between the 
altimeter and the magnetometer indicates that it is possible 
to remove a portion of the magnetic noise due to ocean surface 


gravity waves. 


information be analyzed in the real time domain. Thus, noise 


v 


cancelling must be done real time. 

Spectral analysis or convolution filtering is not practical 
due to thetsignal characteristics of the target anomaly and 
the long length of records required to compute the FFT which 
describes the low frequency signals. 

The magnetic anomaly is a single transient signal imposed 
on the total magnetic field thus making spectral analysis 
impossible. Furthermore, it is possible that the aircraft 
May change direction and flight attitude before a sufficiently 
long record is collected from the altimeter and the magneto- 
meter to determine the coefficients of the Fourier Transform 
£0r convolution filtering. 

In order to keep noise removal in the real time domain the 
peinciple of Adaptive Noise Cancelling (Widrow, et.al., 1975) 


Was applied to the 10 Feb. data. 
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The adaptive noise cancelling was done by a microprocessor 
developed at Naval Oceans System Center (NOSC) and already 
ised in frequency analysis at frequencies higher than those 
analyzed for ultra low frequency magnetic information (0-2 Hz). 
The altimeter (reference signal) contained noise that was corre- 
lated with noise in the magnetometer (primary signal). Each 
Signal was digitized and the algorithm used computed the solu- 
tion which removed the correlated noise from the primary by 
adaptively subtracting the reference signal to converge on 
Minimum power in the primary signal. 

Preliminary results indicate that this method may be useful 
for real time noise cancelling in airborne magnetic data (see 
Figure 7 and Figure 8) while the magnetic anomalies not due 
to wave motion should be relatively unaffected. 

Note in Figure 7 the relative signal levels before and 
after noise cancelling. The signal appears to be cut in gain 
approximately four to one after noise cancelling. However, 
Spectral analysis (Figure 8) indicates that the only power 
reduction was in the frequency band that has been indentified 


as the ocean surface gravity wave noise frequency band. 
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Aircraft magnetic spectrum (A) compared to "noise 


cancelled" magnetic spectrum (B). When altimeter 
was used as noise reference. Each spectrum has 
16 degrees of freedom. Note the band limited 


noise reduction indicated by the cross hatched 
area. 
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VETS ep oUMMARY 


It has been shown that ocean surface gravity waves con- 
tribute significantly to the environmental background noise 
observed in airborne magnetometry. The noise is the result 
of the magnetic field associated with electrical eddy currents 
induced by ocean wave propagation through the earth's magnetic 
field. 

The APN-194 altimeter aboard the P-3C ORION aircraft is 
capable of describing the wave climate below the aircraft sim- 
ultaneously with total field magnetic measurements observed 
by the AN/ASQ-81(V) magnetometer, also aboard the P-3C. Thus, 
the altimeter can be used to infer the wave noise measured by 
the magnetometer. 

Analysis demonstrated that it 1s possible to remove a sig- 
Mmificant portion of the magnetic noise due to ocean gravity 
waves by using the altimeter as a noise reference and sub- 
tracting that processed reference in real time from the mag- 
Netometer time series output. 

Further research with real submarine targets is underway 
in order that the ASW community will benefit from applications 
of this work. Development of a greatly improved magnetometry 


system is possible in the near future. 
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APPENDIX A 
MAGNETOMETER AND ALTIMETER SYSTEM CHARACTERISTICS 


me MAGNETOMETER SYSTEM PARAMETERS - 
AN/ASW-81(V) MAGNETOMETER GENERAL CHARACTERISTICS 


Performance-Magnetometer Group 


Range 

Noise level 

Figure of merit (FOM) 

MAD bandpass (selectable) 


High Pass 


Low Pass 


Gain selection 
(full scale) 


Altitude compensation 


Environmental Limits 


23,000 to 75,000 gammas 
0.016 gamma (peak-to-peak) 
0.24 gamma 


Corner frequencies 0. 
and O..1. Hz 
Rejection slope 36 


Corner frequencies 0.2, 0.4, 0.6, 
and: 2.0: Hz 

Rejection slope 24 dB/octave 

Qed Ow 2g cats. Ly 255 be Oe 20 5 


and 40 gammas 


Barometric type 
Reference altitude, 
Linear range, + 200 


This equipment meets the requirements of military specifica- 
tion MIL-E-5400 within the limits described below. 
Temperature 
(1) Equipment Operating: =54°C to 71°C except the 


Magnetic Detector DT-323/ASQ-81(V) or Towed Body 


TB-623/ASQ-81(V) 
temperature. 


C2) 


which have +55°C maximum 


Equipment Nonoperating: E620 to: 85°C except the 


Magnetic Detector DT-323/ASQ-81(V) or Towed Body 


TB-623/ASQ-81(V) 
temperature. 


Altitude 


Humidity 


which have +71°C maximum 


Sea level to 30,000 feet 


0 to 100 percent 
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Shack 


Vibration 


Input Power Requirements 


Magnetometer Group 


Voltage limits 
Frequency limits 

Power 

Panel Lamp Power 

Voltage limit 

Frequency 

Power 
mMangetometer Pick off Points 

Sienal pack off points in 


pins C-14 and C-15. 


Up to iS ¢ with duration of 11 
+1 milliseconds 


Cunve DV 10" ss “to S00 Hz) except 
Detecting Set Control C-6983/ASQ-81(V) 
or C-9086/ASQ=-81(V) which are curve 


DED Yes. Ge. .9 FO “S00 Ck), 


115/200 volts ac 3-phase power 


130'8> te; 8 voles 
380 to 420 Hz 


150 VA maximum 


0 to 28.5 volts 


de or ac 


10 watts maximum 


the P-3C were Rack C-3, TB-433, 
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mt. ALTIMETER SYSTEM PARAMETERS <= 
AN/APN-194 RADAR ALTIMETER 


Manufactured by Honeywell, Inc., Government and Ared Products 
Division 


Specifications 
tPiGe. accuracy to. 40,000 ft 


cing 


beam width to 3dB point 
-l sec RC constant 
25-35 nano second pulse 
0-25 nsec gate 

Altimeter Pick off Points 


Signal pick off in the P-3C were FS 1005, TB-416, pins 


A-5 and A-6. 
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APPENDIX B 
SPECTRAL WAVE CLIMATE 
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Table B-l1 NOAA Buoy (6003), EB-17 Reported Data for 212 
Jebebe 97 8. 
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APPENDIX C 
DERIVATION OF SURFACE WAVE FREQUENCY SHIFT DUE 
TO MOTION OF AIRCRAFT RELATIVE TO THE WAVE 


The wave equation due to motion of the aircraft is, 


iwt-mx~ iw” t-mx 
e = << 


Where w is the angular frequency of the wave, m is the 
wave number, t is the time and x is the distance from a re- 


aly 


frence point in the direction of wave propagation. The primed 


terms are the values due to motion, 


, 


x “S xP Cucost)< 


Where u is the aircraft velocity andwis the angle clock- 
wise from the horizontal wave motion vector to the horizontal 
aircraft motion vector. 


Then, 


a tot-m(x+(Cucosy)t) 7 2 iw “t-mx 


which leads to 


eivt-mu(cosyp)t ee ae i 


then solving for w” leads to 
@* = W-mu cosy 
when ® = Te m = an and L = & 2s 
equation 26 results 


> =u Lb = 2nfu/ cosy | (26) 
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APPENDIX D 
SPECTRAL VALUES OF AIRCRAFT DATA AND MODELED NOISE 


MAGNETIC CONSTANTS AT 52N, 156N FOR FEBRUARY 1978 


Earth's total magnetic field, Be.= 505807 nt 
° 


Dip angle, I 20 


ES 


THE MODELED MAGNETIC NOISE PRODUCED FROM 
THE “PREDICTED SOWM SPECTRA DATA LEG 1 


FREQUENCY (Hz) AMPLITUDE aT’ paz) 
20? .0002 
03 .0002 
ale) £0024 
14 .0037 
Bails .0029 
pag: /0041 
42/0 .0014 
22 0001 
ve .00064 
6 .0015 
P27 .0028 
29 .0102 
32 .0039 
36 /OO4L 
a37 .0076 
slic .0128 
AAS 0048 
47 .0054 
£56 .0052 
o .0016 
591: .0018 

ioe .00002 


Table D-l.a The Modeled Magnetic Noise Produced From 
the Predicted SOWM Data, Aircraft Course 
PEO Ie 
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2 
PREQUENCY (Hz) AMPLITUDE Ge /Hz) 


eullvall .00087 
«13 s00L2 
14 «0012 
~16 -0037 
os -00LS 
aed ~0014 
«29 .0013 
eicule -0057 
a2 .0088 
ood ~00004 
2912 .0087 
+36 .0069 
Gol -0063 
42 ~0039 
43 .0096 
247 -0034 
pene: -0012 
Pesic) -0022 
68 .0002 
=H sOMOS 7 
od .00019 
29 sO00LS 
i 3 .000018 


Table D-1.b The Modeled Magnetic Noise Produced From 
the Predicted SOWM Data, Aircarft Course 
GCol. 
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Table D-2 FFT Numerical Output of Aircraft Data (Course 270°). 
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